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Abstract. We present results for a multichannel tunneling model that describes point-contact 
spectra between a metallic tip and a superconducting heavy-fermion system. We calculate 
tunneling spectra both in the normal and superconducting state. In point-contact and scanning 
tunneling spectroscopy many heavy-fermion materials, like CeCoIns, exhibit an asymmetric 
differential conductance, dl /dV, combined with a strongly suppressed Andreev reflection signal 
in the superconducting state. For Andreev reflection to occur a junction has to be in the highly 
transparent limit. Here we focus on the opposite limit, namely that of low transparency leading 
to BCS-like dl/dV curves. We discuss the consequences of a multichannel tunneling model for 
CeCoIns assuming itinerant electron bands and localized / electrons. 



1. Introduction 

Point contact spectroscopy (PCS) and scanning tunneling spectroscopy (STS) have been widely 
used to study the electronic properties of heavy-fermion superconductors (HFS). A general 
problem in this field has been the interpretation of tunneling data, which show asymmetric 
conductances in the normal state and significant deviations from the standard BTK formalism 
(Blonder-Tinkham-Klapwijk) [lj. In the past, the BTK formalism has been very successful in 
describing tunneling conductances between metal tips and conventional superconductors, while 
resulting in unphysical parameterizations of the HFS tunneling conductances. 

Here we present results for a multichannel tunneling model between a metallic tip and a heavy- 
fermion superconductor. These results are discussed with respect to the anomalous properties 
observed in CeCoIns [21 El IH El E] j but are readily applied to other heavy fermions (HF) like 
CeCu 2 Si 2 , URu 2 Si 2 and many others El U HD] ■ 

2. Tunneling model 

We model the HF material by two itinerant bands and additional localized surface states, which 
may be caused by broken /-electron bonds at the surface due to the broken translation symmetry, 



The heavy-fermion Hamiltonian %hf represents two bands of itinerant conduction electrons 
with band index a E {light, heavy} and localized electrons near the surface with site index i. 
The operators c' u a (c a .k a ) create (destroy) an itinerant electron with momentum k and spin a 

in band a, while operators fj a (/j CT ) create (destroy) an / electron at site i with spin a. £ a (k) 
are the respective electronic dispersions and Eq is the energy level of the localized / electrons. 
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Figure 1. (Color online) A cartoon of the tunneling processes from the tip of the point contact to 
the heavy-fermion material (localized and itinerant electrons), which are necessary to account for 
the measured asymmetry in point-contact junction conductances and reduced Andreev reflection 
signals. In (a) the localized state appears as a broad resonance at the interface while in (b) the 
localized state forms a sharp surface state, which acts as a resonant tunneling center. 



A simple description of a tunneling experiment is comprised of Hamiltonians for the heavy- 
fermion material, the counter electrode, and the transfer or tunneling processes between them: 
% = %hf + ^electrode + %T- The counter electrode is given by normal conduction electrons 

"^electrode = ^ ] £e O^^k.rr^ka i 
k.cr 



(2) 



and the tunneling Hamiltonian describes all possible transfers 
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In addition to the standard overlap integrals t a between the conduction band in the point contact 
and itinerant heavy-fermion bands there is a finite overlap, ti oc , from the point contact to the 
localized states in the HF. We also account for weak coupling between the localized surface 
electrons and itinerant electrons through scattering terms v a (see Fig. [I]). In general, to get a 
Fano resonance in the conductance one needs interference between different tunneling paths |llj . 
The resulting differential conductance calculated from this model T~Lt will have an asymmetric 
Fano line shape. Figure [T] shows the processes that are active in tunneling between a metallic 
point contact and the HF material. Here we extend the picture of conduction through individual 
quantum channels to a tunneling model to account for point contacts on a HF material. When 
deriving the general tunneling expression, we consider strong overlap between electron states in 
the contact and the HF compound and thus go beyond the strict tunneling limit. 

We calculate the tunneling current through a quantum channel by employing the standard 
non-equilibrium Green's function technique |12 |ll3l[TH ll5j. To further simplify our calculations, 
we make several assumptions: (1) The itinerant microscopic Green's functions are described by 
quasiclassical Green's functions near the Fermi energy. (2) It is essential to keep the full energy 
dependence of the localized Green's function. We assume a single localized level at energy 
e = Eq. (3) Only the heavy electrons undergo a superconducting transition at T = T c , while 
the light electrons remain uncondensed. For details of this model and formalism see Ref. [16J. 



3. Tunneling conductance 

For a multiband tunneling model the differential conductance of a single quantum channel was 
derived in Ref. [IB]. It leads to a generalized Fano expression for the conductance 
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In Eq. (|4j) T> is the transparency of the junction, Eq is the tunneling-renormalized position 
of the localized energy relative to the Fermi level, V is the half-width of the resonance, and 
q = qF + iqB with = B/T. The conventional Fano parameter qp controls the resonance 
shape. The additional parameter qg is present for multiband models only, when tunneling 
through a resonant localized state couples differently to the HF conduction bands (see below). 
The term B adds a Lorentzian to the conventional Fano resonance. 

For notational convenience, we introduce the following parameterization of the microscopic 
parameters t {Khloc] and v {Kl} 



hoc = V%cosa, (5) 



t h = t sin a^/rjh cos 9 t v h = v^/rfa rjh cos 9 V 
t[ = t sin a y/rji sin 9 t ' vi = v^/ijo rji sin 9 V 

where t = t \/N c Nhf, v = vMhf are the effective tunneling elements with the density of states 
at the Fermi level Nshf,c\ 1X1 the heavy fermion (contact) material. The factor r/o is the fraction 
of localized surface states and r]{h,i} give the relative fraction of heavy and light electrons. 
The angles 9t and 9 V give the relative overlap integrals between the direct tunneling and the 
hybridization matrix elements. For simplicity, we assume 9t = 9 V = 9, which results in B = |16| . 
Finally, the angle a quantifies the relative proportion of tunneling into a localized state relative 
to direct tunneling. The four non-zero phenomenological model parameters introduced in Eq. Q 
now depend on the microscopic parameters (t,v, Eo,a,9) and the three fractions ??{o,ft.,i> • 

3.1. Tunneling limit in the normal state 

In the tunneling limit, P < 1, we keep the coupling term v between localized and itinerant 
states at arbitrary strength and obtain 

V « irjiP sin 2 a, (6) 
E w Eq - r] r]ivP sin 2a, (7) 
r w rjolrji v 2 + ^(cos 2 a — rjfv 2 sin 2 a)], (8) 



p , 

qF ~ [1 k(cos z a + r]'fv z sin z a)}. (9) 
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Here r\i = rjh cos 2 9 + rji sin 2 9 is the itinerant fraction. It quantifies the relative weight of direct 
tunneling into light and heavy bands. The parameter v is the effective hybridization matrix 
element at the surface, which is temperature independent. 

In this limit the point-contact spectra measure a all/alV curve with a Fano line shape and 
allow us to directly relate the model parameters to the bulk state of the HF material. At 
temperatures low compared to the coherence temperature of the Kondo lattice, T <S T C ohi 

the 

HF state can be described as a renormalized Fermi liquid [T71 [18] . The density of state factors 
can be related to the ratio of the bare mass m with the effective mass m* as rjh ~ (1 — m/m*) 
and rji ~ m/m* |17j . The temperature dependence of m/m* can be obtained, for example, from 
experiment. For qualitative purposes, we model it as m/m* ~ (T/T co /j) p , with p > 0. The 
factor tjq describes the fraction of localized states and its physical meaning and value are under 
debate, see for example the discussion by Yang [191 12Q] ■ 

We can gain some physical insight by considering three special scenarios on how the measured 
peak width T of the Fano resonance should depend on temperature if it was entirely due to the 
effective mass. For simplicity, we keep only leading order tunneling terms and drop terms 
proportional to t 2 : 

(1) The localized states or moments at the interface are not connected to the HF physics of 
the bulk. In this case, the fraction 770, giving the number of surface states, is temperature 



independent. We obtain 



~2 

T(T) « ^ (1 + 0082(9(1-2—)), (10) 
2 m* 
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where T(T) decreases with increasing temperature, while qF(T) increases, which is inconsistent 
with PCS measurements for CeCoIns, see Ref. |16j . 

(2) Localized states or moments at the interface originating from the Ce 4f states at energy Eq, 
which are the same as in the bulk of CeCoIn§. In this case, it is reasonable to assume that r]o 
is proportional to the number of unscreened moments, i.e., t)q(T) ~ 770 m/m* . We obtain 
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where both T(T) and (If(T) increase with T for T <C T co ^. The increase of T(T) is consistent 
with the PCS measurements, while that of qF(T) is not [16J. 

(3) The dominant temperature behavior ofT(T) is caused by inelastic scattering processes like 
spin fluctuations. Resistivity experiments on CeCoIns show that p(T) may be described by a 
self-consistent spin-fluctuation theory |21j . Similarly, spin- lattice relaxation rates \jT\ ~ T 1 ' 3 
suggest the importance of spin fluctuations [22, 23, 24 . In a spin- fluctuation scenario the 
inelastic broadening increases as T(T) ~ T p , where p > 1, whereas to leading order qp(T) is 
temperature independent. Therefore, broadening due to electrons scattering off spin fluctuations 
is a reasonable interpretation for the temperature dependence of T(T) found in CeCoIns [16J 
and more recently in URu2Si2 HO]. We conclude that the self-consistent incorporation of 
inelastic scattering processes into our calculations, similar to Ref. |25j . is an important next step 
to account for the observed line broadening. 



3.2. Tunneling limit in the superconducting state 

It is necessary to study tunneling in the superconducting state, where Andreev reflection 
processes modify the conductance, because one can extract additional information about the 
microscopic origin of the Fano conductance. The unknown relative ratio of tunneling into heavy 
or light bands can be determined as the dl/dV depends in a non-linear way on 6. At very 
low temperatures, T < T c <C T co h, one can assume that m/m* has saturated to its low-T limit 
~ 1 — rjHF- Thus we write for n defined in the previous section rji ~ t\hf cos 2 9 + ( 1 — tjhf) sin 2 6 
and study the superconducting conductance as a function of 8 and fixing the numerical value 
of i]hf to 0.9. For simplicity, we assume that only heavy electrons become superconducting, 
while light electrons remain uncondensed. In our model, we consider both an s-wave and d- 
wave superconducting order parameter. For the d-wave case we show results for two principal 
orientations. The orientation of the HF crystal lattice relative to the interface normal is either 
antinodal (0°) or nodal (45°) leading to strikingly different tunneling conductances important 
for order parameter spectroscopy of the nodes. 

In Fig. 2 we plot the conductances in the superconducting state calculated in the tunneling 
limit by setting T> = 0.005. The Fano parameters are chosen to fit experimentally observed 
dl/dV curves just above T c = 2.3 K in CeCoIns. Hence we set Eq = 2.01 meV, T = 13.0 meV, 
and qF = —2.16 in agreement with PCS measurements with a Au-tip on CeCoIns, see Fig. 3 
in Ref. |16j . Since the dl/dV curve is only weakly dependent on the tunneling angle a 
in the superconducting state, we set in our calculations a = 7r/20, i.e., stronger tunneling 
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Figure 2. (Color online) Conductances for an HFS in the tunneling limit. The parameters of the 
curves are V = 0.005, E = 2.01meV,r = 13.0 meV,^ = -2.16, a = vr/20, and T = 0.41 K with 
T c = 2.3 K. In panels (a)-(c) we vary the heavy-to-light tunneling parameter 8. The conductances 
are calculated both for an s-wave and d-wave superconducting state, setting A = 0.6 meV. As 
shown, the superconducting dl/dV is very sensitive to the co-tunneling into two bands. The 
inset in panel (a) shows the Fano shaped conductance in the normal state at T > T c with a 
superimposed conductance curve computed at T < T c . 



through the localized state compared to the itinerant bands. In the case of the ideal one-band 
superconducting tunneling model, i.e. 8 = 0, we recover the known dl/dV for all three cases 
considered. In Fig. 2(a) we see the usual BCS superconducting density of states on top of a Fano 
background. In the inset of panel (a) we show the voltage dependence of the dl/dV over a large 
voltage range. In panels (b) and (c) the case of 8 = reproduces the F-shaped dl/dV for the 
0°-junction (panel b), while for the 45°-junction we see the hallmark zero-bias peak (panel c). 
Note that for tunneling into the heavy band only, 8 = 0, the condensate suppresses the normal 
state dl/dV by 100% for voltages below the gap, —A < eV < A. 

Incorporating the additional possibility of tunneling into a band of uncondensed light 
electrons has two main effects on the dl/dV characteristics. First, we see that the coherence 
peaks in the dl/dV curves are reduced from a 100% effect for small 8 to a ~10%-effect, when 
8 > 7r/4. Second, the dl/dV curves are qualitatively changed in shape from their ideal one-band 
8 = appearance. In panel (a) we see that the sub-gap conductance below A can even be 
increased above its normal-state value, resembling a high-transmission sub-gap conductance but 
with a strongly suppressed Andreev reflection signal. Next, turning to the d-wave conductances 
we see that in general, for a fixed set of microscopic parameters, there is a strong dependence 
on the junction orientation. This property may be used as a smoking gun to identify a d-wave 
symmetry. However, if we assume that PCS junctions are made on surfaces with different crystal 
orientations and that these PCS junctions have very different characteristics with different sets 
of microscopic parameters, then examining the curves in panels (b) and (c) may not be as 
distinctive as expected. In this situation one can find cases where the zero-bias peak is strongly 
suppressed making the distinction between principal nodal/antinodal orientations less striking. 
In fact, it may even be difficult to convincingly discriminate between an s-wave and d-wave 
superconducting state based on a few tunneling curves alone. 



4. Conclusions 

We presented results for a multichannel tunneling model for a superconducting heavy-fermion 
material. On the large voltage scale the derived conductances have the features of Fano-like 
dl/dV characteristics. This allows us to extract the microscopic model parameters that describe 
the relevant tunneling processes. Assuming a modified Fermi-liquid state at low temperatures 
and a tunneling contact in the tunneling limit, i.e., low transparency, we find for the normal state 
of CeCoIns that the extracted Fano parameters are most likely describing localized surface states 
and not directly probing the formation of the heavy-fermion state and the temperature behavior 
of the effective mass. In the superconducting state, the calculated dl/dV curves demonstrate 
that additional microscopic information can be obtained about tunneling into heavy vs. light 
bands, which is not possible from studying the normal state alone. Finally, we investigated 
how the dl/dV characteristics depend on the pairing symmetry of the superconducting state. 
We found that in order to reliably extract information from tunneling experiments about the 
symmetry of the superconducting order parameter detailed modeling of the tunneling processes 
is required, which goes beyond a conventional BTK analysis. 
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